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Abstract 

The time reversibility characteristic of Hamiltonian mechanics has long been extended to non- 
Hamiltonian dynamical systems modeling nonequilibrium steady states with feedback-based ther¬ 
mostats and ergostats. Typical solutions are multifractal attractor-repellor phase-space pairs with 
reversed momenta and unchanged coordinates, {q,p) i—> {q,—p) ■ Weak control of the tem¬ 
perature, oc p^ and its fluctuation, resulting in ergodicity, has recently been achieved in a three- 
dimensional time-reversible model of a heat-conducting harmonic oscillator. Two-dimensional cross 
sections of such nonequilibrium flows can be generated with time-reversible dissipative maps yield¬ 
ing aesthetically interesting attractor-repellor pairs. We challenge the reader to find and explore 
such time-reversible dissipative maps. This challenge is the 2015 Snook-Prize Problem. 
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I. TIME-REVERSIBLE NONEQUILIBRIUM FLOWS 


The microscopic gist of the macroscopic Second Law is often pictured as a system’s seeking 
out more phase-space states. So it seems a bit odd that the idealized nonequilibrium steady 
states generated by molecular dynamics behave in the opposite wayi. The (g,p) states from 
such atomistic simulations soon become constrained to an ever-shrinking strange attractor. 
The compensating additional heat-reservoir states are never seen explicitly. Instead they 
are modelled by time-reversible friction coefficients. These coefficients control temperature 
or energy and induce steady-state behavior in the system under study. 

Unlike real life, the underlying Laws of physics are mostly time-reversible, which means it 
is puzzling that physics gives a good accounting of real-world observations. Time-reversible 
mechanical models ( molecular dynamics ) provide clear examples of this conundrum. To 
illustrate the irreversible behavior of such time-reversible models let us consider the simplest 
case, a one-dimensional harmonic oscillator (with unit mass and force constant ) exposed to 
a temperature gradient^>^, T = T{q) . The oscillator’s motion is subject to a time-dependent 
friction coefficient ({p) imposing weak control over the oscillator’s kinetic energy ( K ) and 
its fluctuation, ( ) — { K )‘^ . This oscillator model^ generates a three-dimensional “flow” 

satisfying the three ordinary differential equations of motion : 

q = p; P=-q-C[Ap + Bip^T) ] ; C = ^[ (pVT) - 1 ] + B[ {p^/T^) - ^{p^/T) ] . 

The superior dots in the motion equations indicate comoving time derivatives. 

At thermal equilibrium where T is constant the parameters {A, B) are chosen to promote 
the control of the second and fourth velocity moments^ : 

{A,B) = (0.05,0.32) ^ ( (pVT),(pVt2) ) = (1,3) . 

From the visual standpoint this parameter choice appears to provide an ergodic coverage 
of the oscillator phase space. Liouville’s continuity equation in phase space can be used to 
show that these motion equations are consistent with Gibbs’ canonical distribution ; 

Away from equilibrium the heat-reservoir temperature field T = l + etanh(g) is an explicit 
function of the oscillator coordinate q and e is the maximum value of the temperature 
gradient. 
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FIG. 1: These cross sections can be viewed as generated by time-reversible “dynamical maps”, ad¬ 
vancing the dynamics from one penetration of the C = 0 plane to the next. These ((?,p) phase-space 
cross sections for the conducting oscillator correspond to four values of the maximum temperature 
gradient : e = 0.0, 0.1, 0.2, and 0.3 . Both q and p range from —4 to +A in all the cross sections. 

See Figure 1 for sample numerical (g,p) cross sections of the three-dimensional flow 
in (g,p, C) space. The (g,p) points are plotted whenever changes sign. Notice that the 
equations of motion for this heat-conduction model are time-reversible. To see this start 
out with a solution of the equations forward in time, q{t),p{t), ({t) . Then change the signs 
of the momentum p, (, as well as the direction of time, ( +d/dt —)■ —d/dt ) . These steps 
provide a new “reversed” solution of exactly the same motion equations. 

The nonequilibrium (g,p) cross-sections with e > 0 shown in Figure 1 look like inho- 
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mogenous ( or “ multifractal” ) strange attractors. And they are. In these flows, as in a wide 
variety of time-reversible nonequilibrium steady states, the fractal strange attractors satisfy 
the Second Law of Thermodynamics, with an overall hot-to-cold heat current. The reversed 
flows, topologically similar fractals, are mechanically unstable, with exponentially growing 
phase volume, and with a heat-flow direction violating the second law. These repellor states 
are unobservable numerically due to this ( Lyapunov ) instability, and can only be collected 
by storing, and then reversing, a time series of attractor states^. 

II. EQUIVALENT TIME-REVERSIBLE MAPS^.^^ 

The construction of the flow cross-sections is equivalent to applying a time-reversible 
“map” Ai{q,p) from one penetration of the zero ( plane to the next, {qn+i,Pn+i) = 
M.{qniPn) ■ Such two-dimensional maps provide relatively simple ( two-dimensional rather 
than three ) pictures of the chaos generated by the conducting oscillator. 

Because the penetrations of the differential equations are proportional to the “flux” 
through the sampling plane C = 0 at (g,^, 0) , rather than just the “density”, there are 
unoccupied white “nullclines” in the cross sections wherever the flux C{q,p)f{q,p,0) van¬ 
ishes. 

“Maps” have long been used in dynamical systems theory to illustrate and explain er- 
godicity ( closely approaching all of a system’s states ) and Lyapunov instability ( the 
exponentially-fast growth of small perturbations ). Arnold’s Cat Map>^ and the Baker Map 
are the best-known examples. For clarity let us recall the action of the Baker Map. 

III. THE ROTATED TIME-REVERSIBLE BAKER MAP 

In the dynamical systems literature the discontinuous “Baker Map”, composed of area¬ 
preserving cutting and kneading operations, provides a simple illustration of “chaos”, the 
exponentially diverging growth induced by the mapping of small perturbations. The familiar 
“Cat Map” is another well-known discontinuous mapping of the unit square into itself. 
Both these map types illustrate Lyapunov instability, Sn+i oc 6n , where the proportionality 
constant exceeds unity, leading to divergence. 

Many-to-one mappings, such as Xn+i = axn ; a < 1 , in the presence of chaos are 
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FIG. 2: One iteration of the Baker Map is illustrated here for a grid of 1442 {q,p) points. The 
square shown at the left is first squeezed twofold in the q = p direction and stretched twofold in 
the perpendicular direction. After a cut along the q = p line the two halves are joined together in 
the new arrangement shown at the right. The changed spacings of the points reflect the squeezing 
and stretching induced by the map. 

enough to produce the strange attractors associated with nonequilibrium steady states. 
Even one-to-one volume-preserving mappings, as in the simplest Baker Map above, are 
enough to generate Lyapunov unstable chaos. In the equilibrium Baker Map each half of a 
2x2 square is mapped from a 1 x 2 rectangle to two (1/2) x 2 disjoint rectangles. In order 
to follow our criterion for time reversibility - detailed in the next Section - we use a rotated 
version of the Baker Map, illustrated in Figure 2. The mapping is as follows : 
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-1.5 < q < +1.5 


FIG. 3: The uniform distribution of 100,000 iterates of the Baker Map with an initial point 
{QjP) = (0.5, 0.0) are plotted here. The uniform area conservation of the Map is analogous to the 
equilibrium phase-volume conservation described by Liouville’s Theorem. 

if(q.lt.p) then 

qnew = +1.25d00*q - 0.75d00*p + dsqrt(1.125d00) 

pnew = -0.75d00*q + 1.25d00*p - dsqrt(0.125d00) 

endif 

if(q.gt.p) then 

qnew = +1.25d00*q - 0.75d00*p - dsqrt(1.125d00) 

pnew = -0.75d00*q + 1.25d00*p + dsqrt(0.125d00) 

endif 

The result of 100,000 iterations of this map is a near uniform covering of the 2x2 square, 
as is shown in Figure 3. 
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The dissipative Baker Map^>^, with changes in area as well as shape, is a better model 
for nonequilihrium molecnlar dynamics and statistical mechanics. In those disciplines the 
phase-space density responds to heat transfer. In the oscillator example above the action of 
a heat reservoir with temperatnre T{q) is represented by the time-reversible frictional force, 
—C[ Ap + B{p^/T) ] . The resnlting change in comoving phase volnme 0 follows from the 
continnity eqnation : 

( 0 / 0 ) = {dq/dq) + {dp/dp) + {d(/dC) = -([ A + 3B{p‘^/T) ] . 

The dissipative Baker Map is deterministic, time-reversible, and likewise prodnces mirror- 
image attractor-repellor mnltifractal pairs. These are the same qnalities associated with 
noneqnilibrinm molecnlar dynamics algorithms ever since the early 1970s, thongh they 
passed unrecognized until the 1980s.- 

IV. NEW MAPS - IAN SNOOK PRIZE FOR 2015 

The conducting oscillator is just one example of the deterministic, time-reversible flows 
that represent a nonequilihrium steady state with a chaotic multifractal attractor. The time- 
reversed state, with the momentum, friction coefficient, and the time all changed in sign, 
is an exactly similar mirror-image multifractal structure, an unstable zero-measure repellor. 
More complicated maps with these same properties can be constructed by concatenating 
time-symmetric combinations of reversible mappings { } like 

AllAl2-^3-^4-^3-^2-^l • 

where each of the mappings satisfies the time-reversibility criterion : 

(g,p) = TMTM{q,p) , 

with the four operations performed from right to left. Here T indicates the time-reversal 
operation (g, p) —> (g, —p) . This reversibility criterion states that the sequence of four steps 

- [ 1 ] iterate forward; [ 2 ] reverse velocities; [ 3 ] iterate backward; [ 4 ] reverse velocities 

- returns to the original (g,p) state. The simplest time-symmetric mappings satisfying this 
criterion are shears and reflections.- It is interesting to note that numerical implementations 
almost never return exactly to their initial state after the four-step sequence above. In fact 
the lack of an exact return is a ( somewhat misleading measure of algorithmic accuracy^. 
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The similarities between small-system dynamics and macroscopic dissipative behavior 
motivate the study of relatively simple flows and maps capable of generating the complexity 
associated with irreversible strange attractors. At the same time this complexity often 
exhibits a compelling beauty. The prototypical Cat and Baker maps are relatively simple, 
but their discontinuities are not at all representative of the attractor types shown in Figure 
1 . The time is ripe for a fresh look at such problems. 

The Snook Prize problem for 2015 is to formulate and analyze an interesting time- 
reversible ergodic map relevant to statistical mechanics, in two or three dimensions. It 
is desirable that the map generates a multifractal. Entries can be submitted to the authors, 
to the Los Alamos aryiv, or to Computational Methods in Science and Technology. The 
author(s) of the most interesting entry received prior to 1 January 2016 will receive the 
Prize’s cash award of $500 US . The publisher of Computational Methods in Science and 
Technology ( cmst.eu ) has generously offered an Additional Prize award of $555.55 US 
subject to the Terms and Conditions published in CMST . 
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